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ABSTRACT 

We estimate the fidelity of solar neighborhood (D < 100 pc) catalogs soon to be derived from Pan- 
STARRS astrometric data. We explore two quantities used to measure catalog quality: completeness, 
the fraction of desired sources included in a catalog; and reliability, the fraction of entries corresponding 
to desired sources. We show that the main challenge in identifying nearby objects with Pan-STARRS 
will be reliably distinguishing these objects from distant stars, which are vastly more numerous. We 
explore how joint cuts on proper motion and parallax will impact catalog reliability and completeness. 
Using synthesized astrometry catalogs, we derive optimum parallax and proper motion cuts to build a 
census of the solar neighborhood with the Pan-STARRS 37r Survey. Depending on the Galactic latitude, 
a parallax cut tt/ov > 5 combined with a proper motion cut ranging from /i/c^ >l-8 achieves 99% 
reliability and 60% completeness. 
Subject headings: Data Analysis and Techniques 
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1. INTRODUCTION 

The n ext generatio n of s ynoptic sky surveys (Pan- 
STARRS, EaiseB EQQ3; lsst. iTysonl l2Q02f ) will have the 
necessary astrometric precision and time cadence to mea- 
sure parallax for sources in the solar neighborhood (D < 
100 pc). These parallax measurements will be among the 
most useful signals in building a census of the local stel- 
lar environment. Other popular methods for identifying 
nearby objects (e.g., selection based on proper motion 
or color-magnitude constraints) do not unambiguously 
constrain distance, and require follow-up spectroscopic 
and/or astrometric observations. Furthermore, such se- 
lection criteria are ineffective at discovering sources with 
extreme properties. For example, a bright L dwarf within 
10 pc was only recently discovered in the Sloan Digi- 
tal Sky Survey; it was not found in previous color-based 
searches because of its unusually blue spectral energy 
distri bution (SED) (jSchmidt et all 120101: iBowler et all 
I2010D . Deep, all-sky parallax surveys promise to identify 
nearby objects, regardless of any unusual photometric 
properties, without follow-up telescope investment. 

Currently, the Hipparcos and Tycho catalogs are th e 
only all-sky parallax databases (jPerrvman et al.l 11997). 
While these catalogs are precise, they are also shallow 
( m v 12) and hence biased towards bri ght, high mass 
objec ts. The Pan-STARRS 3ir Survey (Magnier et al. 
2008) is a deep, ground-based synoptic survey that will 
significantly extend the Hipparcos catalogs. The 3n Sur- 
vey will observe 3/4 of the sky 60 times over the course 
of 3 years, and will provide 5-band (grizy) photometry 
for objects in the range 15.3 < m r < 22.© The survey 
has an astrometric error floor of 10 milliarcseconds (mas) 
per exposure, resulting in a stacked parallax sensitivity 
of ~ 2 mas per source. 

This survey will yield a much more complete census of 

beaumont@ifa.hawaii.edu, eugene@ifa.hawaii.edu 
1 Expected detection limits in the gizy filters are 23.4, 22.2, 21.6, 
and 20.1, respectively 



the low mass objects in the solar neighborhood. The first 
Pan-STARRS telescope (PSlQ) saw first light in 2007, and 
the PS1 science mission began on May 1, 2010. 

The 37r Survey will be deeper but less astrometrically 
precise than the Hipparcos catalogs. This poses a unique 
challenge for identifying nearby objects in Pan-STARRS 
data. Distant objects are far more numerous than nearby 
ones and, because of astrometric measurement errors, 
will occasionally display false parallax signals. Thus in- 
terlopers - distant objects with spuriously large parallax 
measurements - will outnumber accurate but small par- 
allax measurements. The focus of this paper is to quan- 
tify this effect, and explore how additional proper motion 
cuts can break the ambiguity between interlopers and le- 
gitimate solar neighbors. 

Interloper confusion is a man ifestation of the Lutz- 
Kelker bias (|Lutz fc Kelkeri[l97l . so we begin with a 
review of that phenomenon in ^JSJ A preliminary anal- 
ysis of the Lutz-Kelker bias suggests that sources more 
distant than 80 pc will be difficult to reliably extract from 
Pan-STARRS data using parallax alone; this is a stricter 
limit than the 100 pc horizon assumed in previous esti- 
mates of Pan-STARR S ' solar neighborhoo d performance 
(jMagnier et all 120081 : iDupuy fc Liul 12009). Thus, in g3] 
we simulate to what extent additional proper-motion 
cuts can suppress interlopers, and what effect these se- 
lection criteria have on catalog completeness. 

We find that cuts based on proper motion and paral- 
lax are more effective than cuts based on parallax alone. 
We estimate that high reliability (99%) catalogs derived 
from Pan-STARRS data will be approximately 60% com- 
plete, and that catalogs with 50% reliability will be 85% 
complete. 

2. THE LUTZ-KELKER BIAS 
ITrumpler fc Weaver! (|1953f ) first noted that parallax 

2 Four telescopes are planned for the long-term Pan-STARRS mis- 
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Fig. 1. — A schematic illustration of the Lutz-Kelker bias. The 
bars show a hypothetical parallax histogram for two classes of ob- 
jects, the more distant of which is also more numerous. When the 
parallax of each source is observed in the presence of noise, the 
resulting histograms spread out into the gaussian curves shown in 
the figure. The group of objects measured to have parallaxes of 
tt' (the dotted line) consists mostly of the more-numerous (red) 
sources. Thus, on average, tt' is an overestimate for sources mea- 
sured to have this parallax. 

measurements tend to systematically underestimate true 
distan ces. This effect was su bsequently analyzed in de- 
tail by iLutz fc Kelkerl (|1973| ) , and has become known as 
the Lutz-Kelker bias (hereafter LK). The LK bias is a 
consequence of the fact that distant stars are more nu- 
merous than nearby ones. To understand this on an in- 
tuitive level, consider the illustration in Figure [1] The 
red and blue bars denote parallax histograms for two 
classes of objects: a large group of distant objects (red 
bar), and a smaller group of nearby objects (blue bar). 
The curves show how these histograms change when the 
parallax is measured in the presence of noise - objects 
scatter away from their true values. Now consider the 
parallax value tt' located midway between each bar (the 
dotted line). There are more objects from the red group 
at this position than from the blue group. Thus, ob- 
jects with measured parallaxes of tt' have, on average, 
actual parallaxes smaller than tt' . This illustration can 
immediately be generalized to continuous distributions of 
stellar distances; when the distribution of stars increases 
with distance, measured parallaxes underestimate true 
distances. 

To quantify this effect, we follow Haywood's (1987) ap- 
proach and use Bayes' theorem to express the probability 
distribution of an object's distance, constrained both by 
a parallax measurement and a-priori knowledge of how 
objects are distributed in space: 

P(D\n) cx P(tt\D)P(D) = N(tt; D-\a K )P(D) (1) 

P(D\tt) is the posterior probability, or the probability 
distribution of an object's true distance, constrained by 
a parallax measurement of tt. P(ir\D) is the likelihood 
of measuring tt, assuming the true distance is D. The 
likelihood is derived from measurement errors, which we 
assume to be Gaussian. In this case, the likelihood is 
given by AT(tt; D -1 , ov), the Normal Distribution with 
mean D' 1 and standard deviation ov, evaluated at tt. 
The prior distribution P(D) is the intrinsic distribution 
of objects as a function of distance, or dN/dD. In an 
observation of solid angle f2, this number is given (to 
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Fig. 2. — The behavior of the Bayesian estimator given by Equa- 
tion [l] assuming a measurement error of cr^ = 2 mas. The three 
curves trace different lines of sight; through the Galactic center 
(solid line), in the Galactic plane along I = 90°, and out of the 
Galactic plane. Top: the mean distance of objects measured to 
have a parallax of 1 /tt; this diverges from the value of 1/tt as that 
number grows. Bottom: the probability that an object's distance 
is >150% the value implied by the parallax measurement. 

within a normalization constant) by 

dN = VLD 2 n(D)dD (2) 

Here, n(D) is the number density of objects as a function 
of D. 

Figure [2] shows the behavior of Equation [TJ assuming 
that av = 2 milliarcseconds (mas); this is the estimated 
maximum precision of the Pan-STARRS 37r Survey (see 
the Appendix). The top panel plots the mean value 
of Equation [T] for three different expressions for n(D). 
They correspond to three lines of sight: one towards the 
Galactic center, one out of the Galactic plane, and one 
away from the Galactic center. Equations for the stellar 
number dens ity are derived from th e Galactic parame- 
ters given in iBochanski et a l. (2010). These three lines 
of sight bracket the behavior of n(D), which increases 
most rapidly towards the Galactic center, and decreases 
most rapidly out of the plane. Note that the mean dis- 
tance of Equation [1] diverges from the parallax estimate 
when 1/-7T > 70 pc. 

The bottom panel plots a similar function: the proba- 
bility that the true distance is more than 150% the value 
implied by the parallax. Again, this probability becomes 
substantial around 1 /tt > 80 pc. Prev i ous estimates 
of Pa n-STARRS yield ([Magnier et al.ll2008t IDupuv fc Liu! 
2009) have assumed that parallax catalogs will be com- 
plete out to 100 pc: a volume twice as large as the 80 pc 
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sphere at which interloper confusion becomes problem- 
atic. 

The limitations of parallax-selected catalogs at moder- 
ate distances raises the question of whether other quan- 
tities might help to increase catalog reliability, without 
excessively compromising completeness. Two promising 
quantities include source color (since the solar neighbor- 
hood is dominated by red, low mass stars) and proper 
motion. In the present work, we estimate the quality 
of solar neighborhood catalogs derived from the Pan- 
STARRS 37r Survey, when both parallax and proper mo- 
tion are used to identify nearby objects. 

3. PROPER MOTION CUTS 

The Bayesian formalism above can, in principle, be 
extended to study how proper motions further constrain 
distance. However, doing so requires defining a prior 
function that realistically models the joint distribution 
of stellar distances and proper motions (that is, P(D) 
is replaced by P(D,fi) in Equation [lj . Furthermore, we 
wish to employ a realistic astrometric error model, which 
depends on source flux. In total, a fully analytic Bayesian 
treatment would require defining a 8-dimensional prior 
distribution (distance, proper motion in two dimensions, 
and flux in 5 Pan-STARRS filters), and marginalizing over 
7 of them. We opt instead of a more straightforward 
numerical approach using simulated Pan-STARRS source 
catalogs. This is more tractable, and allows us to utilize 
existing research into synthesized models of the Galaxy. 

Our approach is as follows. First, we generate catalogs 
of intrinsic source properties (luminosity, color, velocity, 
and distance). We then assign astrometric uncertainties 
and calculate, given a certain set of parallax and proper 
motion selection criteria, the probability that each source 
would be measured to satisfy these cuts. By summing 
these probabilities, we arrive at the expected number of 
sources included in the solar neighborhood catalog. Fi- 
nally, we can compare this number to the actual number 
of nearby objects in the catalogs to measure the perfor- 
mance of the classification strategy. 

Our synthesized catalogs hav e been generated f rom the 
Besancon model of the Galaxy ()Robin et a l. 20030 This 
model accurately represents our current knowledge of 
how stars are distributed in space, velocity, mass, and 
age. Table [T] summarizes the parameters we have used 
when querying the Besancon model service. We have 
generated catalogs for three different positions on the 
sky: (l,b) = (90, 0), (180, 0), and (0, 90). We avoid an- 
alyzing lines of sight near Galactic center for three rea- 
sons: the very high density of sources makes the analysis 
cumbersome; the Pan-STARRS survey will be confusion- 
limited towards this region and hence difficult to analyze 
statistically; and the region is compact on the sky, and 
hence not a significant source of nearby objects. 

For each entry in these catalogs, we calculate the 
(noise-free) parallax, proper motion, and magnitudes in 
the Pan-STARRS (g,r,i,z,y) filter system. We also deter- 
mine the estimated error on each of these quantities. The 
details of this process are covered in the Appendix. 

We next calculate the probability that any object 
would be selected by some cut on parallax and proper 
motion, when those sources are measured in the pres- 

3 http://model.obs-besancon.fr/ 



TABLE 1 

Parameters of the Besanc/on model queries 



Parameter 



Value 



Filter System 

Distance interval 

Small / Large field 

Positions (£, b) 

Solid angle a 

Extinction Law 

Apparent magnitude filtcr b 

Noise 



CFHT-Megacam 
< D < 50kpc 
Small 

(90,0), (180,0), (0,90) 
4, 20, 1000 deg 2 
mag kpc — 1 
r < 26 
none 



a Solid angles were chosen to generate large but manage- 
able catalogs, containing 4-6 million entries, 
k This range is wide enough to include all stars detected 
in the Pan-STARRS filter system. 

ence of realistic noise. We assume that all measurement 
errors are Gaussian, so that the probability of measur- 
ing a parallax greater than some threshold 7r t is obtained 
directly from the Gaussian cumulative distribution func- 
tion. We also consider the probability that the proper 

motion (that is, [i = \J fJ% + A*y) exceeds some threshold 

[it- When jj, x and \i y are measured with uncertainty cr^, 
this probability is given by 



Pfa > ih) = \ 



x, 2, -5- ) dx 



(3) 



where % 2 (x, v, A) is the non-central chi-square distribu- 
tion, with v degrees of fr eedom and no n-centrality pa- 
rameter A, evaluated at x (|Patnaiklll949f ). 

3.1. Defining Completeness and Reliability 

The task at hand is to choose some combination of se- 
lection criteria on w, fi 7 and the errors on these quantities, 
to optimize the solar neighborhood catalog. To do so, we 
define the quantities reliability and completeness to de- 
scribe catalog fidelity. These are application-dependent 
concepts, and different definitions lead to different "opti- 
mal" selection criteria. We wish is to assess the utility of 
Pan-STARRS data in identifying objects within the solar 
neighborhood, so we use the following operational defi- 
nitions for completeness and reliability: 

The completeness of the Pan-STARRS solar neighbor- 
hood catalog is defined to be the fraction of observable 
objects within 100 pc that are included in the catalog. An 
object is said to be observable if it is detected but not 
saturated in at least one filter at each epoch. The value 
of 1 00 pc is chosen for consistency with previous stud- 
ies (jDupuv fc Liul l2009). and represents a ~ 5a parallax 
detection for bright sources (see the Appendix) . The re- 
liability is defined as the fraction of the catalog occupied 
by sources with true distances less than 200 pc. Note the 
"buffer" from 100-200 pc; this is designed so that only 
sources substantially more distant than 100 pc degrade 
the reliability. Including a source within 200 pc does not 
compromise a catalog's census of nearby objects. How- 
ever, a source beyond 200 pc displays a parallax signal 
measured at < 2.5er significance, and hence would not be 
considered a "detection" by most criteria. 

In summary, the completeness gives the probability 
that a nearby source is in the catalog (provided it was 
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detected photometrically); the reliability gives the prob- 
ability that a random catalog object is indeed in the solar 
neighborhood. 

We have subjected the catalogs described in the previ- 
ous section to several filtering strategies, and measured 
the completeness and reliability for each strategy. Each 
strategy combines three different cuts: an absolute par- 
allax cut (1/tt < 200 pc, or n > 5 mas), and a cut in rel- 
ative parallax and proper motion (ir /ov > x, h/<j 11 > y). 
We sampled a grid of strategies, ranging from < x < 7 
in steps of 1 and < y < 30 in steps of 5. We also 
sampled a sub-region 2.5 < x < 5 in steps of 0.2 and 
< y < 9 in steps of 0.5. In each simulation, any object 
meeting all criteria (absolute parallax, relative parallax, 
and relative proper motion) is labeled as a nearby source. 

4. RESULTS 

Figure [3] plots the completeness (greyscale) and reli- 
ability (red contours) for different selection criteria and 
fields. Consider first the completeness for each field. As 
expected, completeness decreases monotonically for in- 
creasing cut strength. Note that the b — 90° field is less 
sensitive to proper motion cuts. This is because proper 
motion measurements for lines of sight at high latitudes 
trace motion parallel to the plane, while measurements 
at low latitudes are more sensitive to vertical motions. 
Because velocities in the plane are faster than those out 
of the plane (both in absolute terms and relative to the 
solar motion), objects towards b = 90° display higher 
proper motions, and proper motion cuts have less of an 
effect on completeness in this field. 

The reliability varies more substantially between the 
three fields. As was seen in Section [2] lines of sight closer 
to the Galactic center display lower reliability for a given 
selection strategy; this is due to increased contamination 
from distant objects. Especially in the range ir/a^ > 
x, 2 < x < 5, more aggressive proper motion cuts in- 
crease reliability. Furthermore for [ijo^ > y, y < 10, 
these cuts have little impact on completeness. 

Figure [3] presents an opportunity to identify optimal 
filter strategies, by maximizing completeness along con- 
tours of constant reliability (or vice- versa). Figure |4] dis- 
plays the results of this optimization. It shows the maxi- 
mum completeness obtainable for any given reliability. 
Also shown in dotted lines are the optimal strategies 
without proper motion cuts. For most values of relia- 
bility, proper motion cuts increase the completeness by 
5 — 10%. Because non-optimal strategies fall below these 
curves - and because the curves are monotonic - Fig- 
ure 2] also displays the maximum reliability for a given 
completeness. Especially for high completeness values, 
stronger proper motion cuts improve reliability by >10%. 

Proper motion provides the greatest benefit to lines of 
sight in the Galactic plane. This is where interloper con- 
fusion is most severe, and where additional constraints 
help to resolve the confusion. Conversely, proper mo- 
tion filtering has little impact on lines-of-sight out of the 
plane. 

Table [5] lists the specific strategies that lead to the 
curves in Figure 2J In general, high reliability strategies 
combine 5a parallax cuts with 1-8 a proper motion cuts, 
depending on Galactic latitude. 

5. DISCUSSION 




1 2 3 4 5 6 7 

7i/<y„ 




1 2 3 4 5 6 7 

71 / a„ 




1 2 3 4 5 6 7 

it /a, 



Fig. 3. — The completeness (greyscale) and reliability (red con- 
tours) for different selection criteria and lines-of-sight. Proper mo- 
tion cuts have the greatest benefit for tt/cf-k > x, /j/o> > V m 
the range 2 < x < 5 and y < 10. In this region, more aggressive 
proper motion cuts improve reliability, and has little impact on 
completeness. 
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Fig. 4. — The completeness and reliability of optimal strategies. 
Non-optimal strategies fall below these curves. The solid lines trace 
strategies using both parallax and proper motion selection, while 
the dotted lines trace parallax cuts only. The slight jaggedness in 
the curves is an artifact of our finite sampling of strategies. 



TABLE 2 
Optimum Filter Strategies 



Reliability 


Max Completeness 






t = 90° b = 0° 


0.50 


0.83 


2.9 


6.2 


0.60 


0.81 


3.2 


5.5 


0.70 


0.79 


3.4 


5.5 


0.80 


0.76 


3.7 


5.5 


0.90 


0.72 


4.1 


1.5 


0.95 


0.68 


4.3 


7.1 


0.99 


0.59 


4.9 


8.0 


i = 180" 6 = 0° 


0.50 


0.83 


2.9 


4.2 


0.60 


0.82 


3.1 


4.0 


0.70 


0.79 


3.4 


4.0 


0.80 


0.77 


3.6 


4.0 


0.90 


0.73 


4.0 


4.0 


0.95 


0.69 


4.3 


1,1 


0.99 


0.59 


4.9 


5.4 


b = 90" 


0.50 


0.86 


2.8 


3.5 


0.60 


0.84 


3.0 


3.0 


0.70 


0.82 


3.3 


3.0 


0.80 


0.80 


3.6 


2.5 


0.90 


0.76 


4.0 


2.0 


0.95 


0.72 


4.3 


2.0 


0.99 


0.63 


5.0 


1.0 



Figure 2] demonstrates that parallax and proper mo- 
tion cuts lead to more complete and reliable solar neigh- 
borhood catalogs than strategies based on parallax cuts 
alone. This is true mainly because proper motion se- 
lection is effective at identifying the false parallax sig- 
nals discussed in Section [5J While parallax provides the 
best quantitative constraint on distance, proper motion 
is effective at discriminating between distant sources and 
those in the solar neighborhood. 



The main disadvantage of proper motion selection is 
that, unlike parallax, there is not a one-to-one relation- 
ship between proper motion and distance. Because of 
this, proper motion selection necessarily introduces bi- 
ases into the extracted catalog. For example, inter- 
lopers in a proper-motion selected catalog will be bi- 
ased towards population II stars, on higher-velocity or- 
bits through the Galactic potential. In addition, nearby 
sources on orbits similar to the sun's will be culled. 

An alternative strategy would incorporate color cuts 
to suppress interloper confusion. This would affect the 
resulting catalog in a different way; color cuts would bias 
against unusually red or blue sources. Depending on the 
science application, this bias might be less desirable than 
the bias introduced by proper motion cuts; extreme col- 
ors imply unusual and potentially interesting physical 
conditions in the (sub) stellar atmosphere. 

Finally, we note that while the above simulations illus- 
trate how to design complete and reliable source catalogs, 
the quantitative performance of each strategy depends 
on the specific error properties of the data. We have 
assumed sensible noise parameters describing the Pan- 
STARRS 3tt Survey, but these parameters are likely to 
evolve as the survey matures. In addition, different scien- 
tific motivations will lead to different operational defini- 
tions for completeness and reliability. Both of these fac- 
tors will change the detailed performance of each strat- 
egy. Thus, we recommend repeating simulations like 
these when deriving the final solar neighborhood cata- 
logs from Pan-STARRS data. 

6. CONCLUSION 

In this report we have investigated to what extent the 
Pan-STARRS 37r Survey will be able to reliably identify 
nearby objects (D <I00 pc). We have argued that iden- 
tifying such objects by parallax alone is not optimal; such 
a strategy is substantially contaminated by interlopers, 
due to the fact that nearby sources are so heavily out- 
numbered by distant ones. Because of these interlopers, 
Pan-STARRS parallax measurements will be ineffective at 
identifying nearby objects more distant than 80 pc (as- 
suming a n = 2 mas; Section [2]). 

Additional proper motion cuts are able to further dis- 
criminate between nearby and distant objects. We have 
derived a set of optimal strategies incorporating both 
parallax and proper motion cuts, which maximize com- 
pleteness for a given level of reliability (and vice- versa). 
We estimate that 99% reliable catalogs will be 60% com- 
plete, while catalogs at 50% reliability will be 85% com- 
plete. 

Compared to other possible cuts (e.g. selection on 
color) we suspect that proper motion cuts are less suscep- 
tible to biasing against finding as-yet unidentified classes 
of objects in the solar neighborhood; these putative ob- 
jects may well have unusual colors, but would mostly 
display typical proper motions. 

We thank T. Dupuy for helpful discussions about Pan- 
STARRS sensitivity calculations, and A. Robin for discus- 
sions about the Besancon model. 



6 



Beaumont & Magnier 



TABLE 3 
Pan-starrs System Parameters a 



Area 
(m 2 ) 


^fwhm 

arcsec 2 


^pixel CT read -^cpoch ^baseline 

arcsec 2 electrons px — 1 years 


1.73 


0.8 


0.26 5 12 b > c 3.5 C 



a Data from Ma gnier et al.H200Sl and lOnaka et alj|2008t 

Epochs per filter. Due to gaps in the focal plane array, the average 
object will be detected in 80% of the epochs. 
c Specific to the 3tt Survey 



TABLE 4 
Pan-starrs Filter-Dependent 
Parameters' 1 



Band t cxp b $ C * s ky d 
(s) (10 9 s" 1 ) (s- 1 ) 



g 
r 
i 

z 

y 



60 8.18 17.3 

38 10.2 56.5 

30 9.60 122 

30 6.54 195 

30 3.20 225 



a Data from IDupuv fc Liull2009l) 

b For the 3tt Survey 

c Photoclcctron rate from a zero magni- 
tude source. 

^ Photoclcctron rate for the sky. As- 
sumes that each source measurement 
samples an effective angular area of 
4-7T(7^ Rf , appropriate for PSF fitting 
(Mighcll 2003). 

APPENDIX 

SYNTHESIZING PAN-STARRS MEASUREMENTS FROM THE BESANCON CATALOGS 

The Besangon model generates realistic synthesized catalogs of sources one might observe at a given location on the 
sky. Because we are interested in the specific performance of the PS1 3n Survey, we determine apparent magnitudes 
in the Pan-STARRS (g,r,i,z,y) system and, from these quantities, the corresponding astrometric errors. Relevant pa- 
rameters desc ribing the Pan-STARRS PS 1 telescope, Gigapix el Camera, and 37r Survey are provided in Tables [3] and H] 
(jOnaka et alJ l2008: M agnier et alJl200l IDupuv fc Lh]f200l . 



The Besangon catalogs return magnitudes in the CFHT Megacam photometric system. The Pan-STARRS and Mega- 
cam g filters are similar, so we accept the g magnitude included in the catalog by default. To determine fluxes in the 
other bands, we derive colors for each star by convolving a re presentative stellar s pectrum with the Pan-STARRS filter 
throughput p rofiles. We follow the same process discussed in IDupuv &: Liu ([2009( 1. and the spectra we use come from 
iPicklesI Jl998). We approximate the white dwarfs in the Besangon catalogs as blackbodies, and AGB stars as class I 
M5 giants. We manually model interstellar reddening via an extinction rate of 0.7 mag Ay kpc -1 , and an extinction 
curve from Cardelli et al. (1989). Lines of sight are only reddened while they are within 300 pc of the Galactic plane, 
since this is where most Galactic dust is confined. Note that the photometric quantities themselves are not directly 
used in our analysis; they are used merely to derive astrometry errors, and to determine detection and saturation 
limits. 

To determine the Poisson noise associated with each source, we convert apparent magnitudes to photo-electron 
counts given the parameters in Tables [3] and |4j We assume photometry is performed via PSF fitting, such that all 
photo-electrons from the source are used when estimating the Poisson source noise. 

We include two additional sources of noise: Poisson noise fr om the sky, an d read noise. For PSF-fitting photometry, 
the effective area of sky sampled for noise estimates is A-Ka^ s{ (|Mighellll2003t l . The total number of sky photoelectrons 
measured with each source is given in Table |4] By adding the source-, sky-, and read- noise terms in quadrature, we 
arrive at the photometric signal-to-noise ratio for each measurement. 

We model two sources of astrometric error. The first is due to the broadened point spread function (PSF) induced 
by the atmospher e, which we t ake to be Gaussian. The ability to determine the center of this Gaussian is given 
approximately by ()Mighelll l2003') 

" = it < A1 > 



where <7 ps f is the width of the PSF and SNR is the signal-to-noise ratio derived above. 
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The second astrometric noise source reflects the f act that the PSF cen ters are randomly displaced from their true 
centers because of turbulent atmospheric refraction (|Christou et al.ll2003f) . Most of this refraction is correlated among 
nearby stars, and can be corrected for during data reduction. However, the differential refraction between nearby 
sources (known as tilt isoplanatism) is not correctable. This adds a random offset to each stellar position, and 
sets an error floor (ct 2 ) for each exposu re; this error domin ates <j\ for bright sources. The Pan-STARRS error floor 
is approximately 10 mas per exposure (M agnier et al.l l2008f ). Assuming that this is also Gaussian distributed and 
statistically uncorrelated with o\ , the net astrometric error per epoch is given by 

Equation IA2I gives the positional error for a single observation, in a single filter. We are mainly interested in the 
parallax and proper motion uncertainty, given several position measurements with uncertainty er. From the Central 
Limit Theorem, for sufficiently large TV we expect 

a a 

<Jk (X —= (Tp oc -j= (A3) 

V A ibaselineV N 



where A is the number of observations, a is given in equation I A21 and ^baseline is the time baseline of the observations 
in years (since proper motion is reported in units of angular displacement per year). As Table [3] indicates, each field 
will be observed 12 times in each of 5 filters, over a period of 3.5 years. Because gaps fill 20% of the focal plane array, 
the average source within this field will be observed 9.6 times. 

We determine the constants of proportionality for Equation IA3I in a Monte Carlo fashion, by fitting parallax and 
proper motion to many collections of synthetic observations. While a n varies somewhat with position on the sky, 
the two constants of proportionality are given approximately by 1.3 for cv and 3.5 for The proper motion error 
corresponds to the uncertainty along each direction. 

To summarize, then, equation A3, combined with the constants of proportionality given above, give the per-filter, 
stacked error in parallax and proper motion. These errors can be further reduced by stacking the signal from each 
filter. We thus compute the standard error of the weighted mean of the five per-filter parallax estimates. This number 
finally gives the parallax and proper motion uncertainty for each object in the catalog. 

As an example, consider a bright source -C 1 for every filter). Then, from equation IA21 a = er 2 = 10 mas, 

per filter, per exposure. Using equation IA31 the net per-filter parallax uncertainty is a v = 1.3 x 10/V9.6 = 4.2 mas. 
Since each filter has the same uncertainty in this case, the standard error on the weighted mean of the five per-filter 
measurements is ov = 4.2/V5 = 1.88 mas. Note that, in general, the five per-filter uncertainties need not be the same 
- SNR often varies across filters. 
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